Polarization Dependence of Raman Spectra in Strained Graphene 
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The polarization dependences of the G, D, and 2D (G') bands in Raman spectra at graphene 
bulk and edge are examined theoretically. The 2D and D bands have different selection rules at 
bulk and edge. At bulk, the 2D band intensity is maximum when the polarization of the scattered 
light is parallel to that of incident light, whereas the D band intensity does not have a polarization 
dependence. At edge, the 2D and D bands exhibit a selection rule similar to that of the G band 
proposed in a previous paper. We suggest that a constraint equation on the axial velocity caused by 
the graphene edge is essential for the dependence of the G band on the crystallographic orientation 
observed in the bulk of strained graphene. This is indicative of that the pseudospin and valleyspin 
in the bulk of graphene can not be completely free from the effect of surrounding edge. The status 
of the experiments on the G and D bands at the graphene edge is mentioned. 
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I. INTRODUCTION 

Since the early stage of the research on graphene, char- 
acterization of a sample has been the central issue and 
Raman spectroscopy has been playing a major role in 
characterizing a sample. 1 For example, the 2D (G') band 
in Raman spectra is useful in distinguishing a monolayer 
from few-layer graphene stacked in the Bernal configu- 
ration, and the appearance of the D band indicates 
that an intervalley elastic scattering of a photo-excited 
electron is activated by defect.— ~— Another advantage of 
Raman spectroscopy, besides the characterization of a 
sample, is that Raman spectra can include detailed in- 
formation on the wave function of the electron. 

Raman process concerns with photon, phonon and 
their mutual interaction through the electrons. Because 
the electron-photon and electron-phonon interactions in 
graphene are relevant to pseudospin and valleyspin^ the 
Raman spectra are capable of retrieving information on 
the pseudospin and valleyspin. An interesting point here 
is that the pseudospin and valleyspin are sensitive to the 
presence of graphene edge<2 As a result, we can have a 
selection rule specific to the graphene edge. For exam- 
ple, it is known that only the armchair edge enhances the 
D band intensity and that the intensity depends on the 
angle between the armchair edge and the polarization of 
incident (scattered) laser light 

In a previous paper, we proposed a selection rule for 
the G bandji&i^ This selection rule states, for exam- 
ple, that the Raman intensity is enhanced when the po- 
larization of Raman laser is parallel (perpendicular) to 
the armchair (zigzag) edge. This prediction has been 
supported by recent experiments of Cong et alr^ and 



Begliarbekov et ai^ Their experiments illustrate that 
the G band intensity exhibits the anomalous polarization 
dependence at graphene edges which is different from the 
polarization dependence at the interior (bulk). Their re- 
sults could be naturally explained in terms of the special 
behavior of the pseudospin and valleyspin near the edges 
of graphene. 

In this paper, we explore selection rules for the D and 
2D (G') bands at bulk and edge. Since the 2D band is 
a prominent peak in Raman spectra, the selection rule 
must be useful in extracting more information on the 
pseudospin and valleyspin from the Raman spectra. In 
addition, we examine the G band in strained graphene 
as the application of the selection rule for the G band. 
It is known that strain splits the G band into two sub- 
bands called G + and G~, and that the Raman intensity 
of each subband has a crystallographic orientation depen- 
dence i 17 ' 18 Our result suggests that the crystallographic 
orientation dependence observed in the bulk of strained 
graphene is relevant to the selection rule of the G band 
for the graphene edge. The pseudospin and valleyspin in 
the bulk of graphene seem to be not completely free from 
the effect of surrounding edge. 

This paper is organized as follows. In Sec. [TT] we de- 
rive the selection rule of the G band for the graphene 
edge in a unified manner. In Sec. IIIII we apply the con- 
straint which is essential for the selection rule of the G 
band to explaining the crystallographic orientation de- 
pendence of the G band Raman intensity observed in 
strained graphene. The selection rules of the D and 2D 
bands are proposed in Sec. IIVI Discussion and summary 
are given in Sec. [V] 
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II. QUICK OVERVIEW OF SELECTION RULE 
FOR G BAND 



In this section we reproduce the selection rule of the 
G band for the graphene edge obtained in a previous pa- 
per^ by employing an approach based on two velocities 
associated with the gauge helds for photon and phonon. 
This new approach can help us to recognize strange simi- 
larity between the zigzag and armchair edges. This simi- 
larity is represented by the condition Eq. (fTTj) or Eq. (fT6| . 
As we will show in Sec. this condition is necessary to 
explain recent experiments showing that the G band ex- 
hibits a polarization dependence on the crystallographic 
orientation of strained graphenei 17 ! 18 



A. Two Velocities 

Let us begin with the Hamiltonian including photon 
field A and phonon fields A q and <j)^ 



and that of v q as 



H 



a ■ (p - eA + A q ) (f>a x 

<p*<j x a' ■ (p - eA - A q ) 



(1) 



where p = — iV is the momentum operator, a a (a — 
0,x,y,z) is the pseudospin, a = (a x ,a y ), and <x' = 
(—cr Xl (J y ). The phonon field <f> (A q ) gives rise to an in- 
tervalley (intravalley) scattering. In Eq. (|T|), we have 
adopted units in which fi = 1 and vp = 1, and omitted 
the position dependence in the variables A q and cf> be- 
cause we are interested in the T and K points phonon 
modes. 

From Eq. |T]), we define two velocity operators, v and 
v q , as follows: 



1 dH 
V = ~~edA 
dH 



dA q 



a 
a' 

■ 
-er' 



(2) 



The operator v is nothing but the usual velocity operator 
that couples to an electro-magnetic gauge field A. Note 
that an electromagnetic current is given by multiplying 
—ev-p and v together, and the unperturbed Hamiltonian 
is written as H — vp. The velocity v q is distinct from v 
for the sign in front of er'. It may be appropriate to call v q 
an axial velocity because v q has some analogy to the axial 
current in quantum electrodynamics. In graphene, the 
axial velocity v q couples with a lattice deformation. Note 
that phonon is one example of a lattice deformation and 
a general lattice deformation such as ripples and edges 
can be represented by an axial gauge field A q (r)£t~— 

By using r a (a = 0, x, y, z) to represent the valleyspin, 
we write the component of the velocity v as 



<JxT z , 



(3) 



Ca;To, 

(TyTg. 



(4) 



Note that these two velocities v and v q are related with 
each other via the pseudospin a z as 



-ieijvj, 



(5) 



where i,j 6 y} and is antisymmetric tensor satis- 
fying e xy = 1, e yx = -1, and e xx = e yy = 0. In the fol- 
lowing subsections, by examining the effect of the zigzag 
and armchair edges on the two velocities v and v q , we 
derive the selection rule for the G band. 
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FIG. 1: (a) The coordinate system represents the crystal ref- 
erence system. The zigzag edge is parallel to the a;-axis, while 
the armchair edge is parallel to the y-axis. (b) A graphene 
system with rectangle shape enclosed by the zigzag and arm- 
chair edges. 



B. Zigzag Edge 

First of all, the electronic velocity v normal to the 
zigzag edge must vanish i 22 i 23 By taking the zigzag edge 
along the x-axis [see Fig.[TJa)] , we thus have the condition 



(Vy) = 0. 



(6) 



Here, (O) denotes the expectation value of the operator 
O with respect to the standing wave near the edge. A 
general wave function can be represented by 



*(r) 



* K (r) 
*K'(r) 



(7) 



where "^(i 1 ) [$K'(r)j is the two-component wave func- 
tion for an electron near the K [K 7 ] point. The two- 
component structure corresponds to the pseudospin. 
Note that the zigzag edge is not the source of intervalley 
scattering ] 14 i 24 As a result, the standing wave near the 
zigzag edge is written as 



*(r) 



* K (r) 






*K'(r), 



(8) 
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Thus, Eq. © leads to 



^ K (r)ay^ K {r)dxdy = 0, 



/ ^ K ,(r)a y ^ K '(r)dxdy = 0. 
Js 



(9) 



Namely, the condition (v y ) — must be satisfied inde- 
pendently for the K and K' points. Because 



(r + r z )*(r) - 2 



(r - r,)*(r) = 2 



*K(r) 




*K'(r) 



(10) 



Eq. © is possible only when (oj,(to±t*)) = [i.e., (v y }± 
(v q ) = 0] is also satisfied. Hence, we get a constraint for 
the axial velocity as 



(t# = 0. 



(11) 



We will show below that the selection rule for the G band 
arises from the two conditions given by Eqs. © and (JTT|) . 

The Hamiltonian contains the electron-phonon (el-ph) 
interaction, A q • v q , and the constraint Eq. ([TO shows 
that A q component does not have a nonzero el-ph matrix 
element. Because the vector A q is pointing perpendicular 
to the optical phonon eigen vector ) 14 ' 25 i 26 A q component 
corresponds to the optical phonon mode whose displace- 
ment vector is parallel to the zigzag edge, u x , which we 
called the longitudinal optical (LO) phonon mode in a 
previous paperj 13 i 14 Thus, Eq. pT|) shows that the LO 
mode is not a Raman active mode near the zigzag edge^ 4 ® 
Only transverse optical (TO) phonon mode A q (u y ) can 
be Raman active. 

Moreover, since the pseudospin a z changes the wave 
function from symmetric (bonding) to anti-symmetric 
(anti-bonding), the optical transition amplitude is pro- 
portional to the expectation value of va z . The optical 
transition does not take place when the polarization of 
the incident (or scattered) laser light is parallel to the 
zigzag edge (or the x-axis) because the corresponding 
optical matrix element vanishes as 



A x (v x <j z ) = -iA x (a y T z ) = -iA x (v^) = 0, 



(12) 



due to Eq. (fTTjl . Similarly, the phonon softening (Kohn 
anomaly) is relevant to the expectation value of v q a z . 
The phonon softening is absent for the unique Raman 
active TO (A q ) mode because 

A*{t£<r,) = -iA%{v y ) = 0, (13) 

due to Eq. ©. Note that the LO mode (A q ) can un- 
dergo a phonon softening effect, however, the LO mode 
is invisible to Raman spectra. 

In summary, the polarization of the laser light should 
be perpendicular to the zigzag edge in order to have a 
Raman intensity, and the corresponding Raman active 
mode is the TO mode which is free from the phonon 
softening effect. It is important to recognize that this 
selection rule is a consequence of the conditions Eqs. © 
and Ijlip for the electronic and axial velocities. 



C. Armchair Edge 

The armchair edge can be examined in a manner sim- 
ilar to that for the zigzag edge. The electronic velocity 
normal to the armchair edge must vanish. By taking the 
armchair edge along the y-axis (see Fig. [Ha)), we have 
the condition 



0. 



(14) 



Note that the armchair edge is not the source of intraval- 
ley scattering and preserves the pseudospin under an in- 
tervalley scattering] 14 i 24 As a result, we obtain [i — x, y) 

^ K (r)ai^ K {r)dxdy = / ^ K ,(r)ai^ K '(r)dxdy. (15) 
Js 

This equation is equivalent to the conditions (a x T z ) = 
and {(J y r z ) — 0. The former condition (<j x t z ) = is 
nothing but Eq. (fl4]) . while the latter one (<J y T z ) = 
corresponds to 



0. 



(16) 



From this condition it is straightforward to show that 
only the LO mode A q (u y ) is Raman active mode at 
the armchair edge. The optical transition does not take 
place when the polarization of the incident laser light is 
perpendicular to the armchair edge because 



A x {v x a z ) = -iA x (v q ) = 0, 



(17) 



due to Eq. (fT6"j) . Furthermore, we see by using Eq. (|T4"|) 
that the Raman inactive TO mode {A q ) does not undergo 
the phonon softening because 



A q {v%a z )=iA y -{v x 



= 0. 



(18) 



The Raman active LO mode (A q ) can exhibit a phonon 
softening effect. To summarize the selection rule of the G 
band for the armchair edge, the polarization of the laser 
light should be parallel to the armchair edge in order to 
have a Raman intensity, and the corresponding Raman 
active mode is the LO mode which undergoes a phonon 
softening effect. 

The zigzag and armchair edges are distinct concern- 
ing the usual velocity v as (v y ) — and (v x ) = 0, re- 
spectively. Note, however, that the zigzag and armchair 
edges are not distinct with respect to the axial velocity 
v q . The constraint for the axial velocity at the armchair 
edge is the same as that at the zigzag edge as shown by 
Eqs. (fTTjl and (fT6")) . although the origins of Eqs. (TTT1) and 
(fT6")) are totally different. For the zigzag edge, Eq. pT|) 
is satisfied both for the K and K' points since the zigzag 
edge is not the source of an intervalley scattering. In 
other words, the zigzag edge affects only the pseudospin 
as (<7 y ) = Of 1 ^ The zigzag edge is irrelevant to the val- 



leyspin. In contrast, Eq. (jT6|) is satisfied because the 
armchair edge is the source of an intervalley scattering 
and preserves the pseudospin. The armchair edge affects 
only the valleyspin and is irrelevant to the pseudospin. 
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III. UNIAXIAL STRAIN 

In this section we apply the selection rule of the G band 
for graphene edge to the G band in strained graphene. It 
is known that strain lifts the degeneracy of the G band, 
so that the G band splits into two subbands called G + 
and G~. The eigenvectors for the atomic displacements 
are perpendicular to the direction of strain for the G + 
mode, and parallel to the G~ mode i 17 ' 18 By defining A q 
(Al) field for the G+ (G~) mode, we see that A^ is 
parallel to the direction of strain, while A_ is perpen- 
dicular to it since A q is pointing perpendicular to the 
corresponding optical phonon eigenvector. Recent ex- 
periments by Huang et altH and Mohiuddin et alr^ for 
bulk of strained graphene show that the Raman inten- 
sity of each G + and G~ mode has the crystallographic 
orientation dependence. We suggest that the crystallo- 
graphic orientation dependence of the G band Raman 
intensity observed at bulk can be attributed to the con- 
straint (v q ) = obtained in the previous section for the 
graphene edge. The difference between the bulk and edge 
is pointed out with respect to the crystallographic orien- 
tation dependence of the G band Raman intensity and 
Kohn anomaly effect. 



(b) 




photo-excited electron is a resonant state. Note that this 
effective interaction is for the bulk where the electron 
propagator does not depend on the pseudospin. Owing 
to Eq. ([3]), Hq can be expressed in terms of the axial 
velocity operator only. Moreover, by using the properties 
of the axial velocity operator, vfvf = gqTq (i = x, y) and 



v q v q 



v q v q , we rewrite %q in the forrn^ 



/ _in _out 

2 \ x x 

(An ^out 
~ \ € x € v ' 



^in^out 



)(Ay« + A«v%). 



in out 
t y e x 



(20) 



Note that Hq is linear in v q . This feature is unique to the 
G band and is not seen in the case of the D and 2D bands 
as we will show later. By introducing the angles, 9 ln , O ut, 
and (fi s as shown in Fig. [2j one has e™ = e ln cos(6>i n + ip s ), 



'sin(6> in + ip s ), e° 
?out +¥>s), and 



3(#out + <p s ), 



A q = Al cos ip s 
A* = Al shnp 3 



- Al sin (p s , 
Al cos ip s . 



(21) 



By putting these into Eq. (|20p . we get 
H G 



{v*A\ - v q Al) cos(0 in + 9 out + <p s ) 



(v q Al+v q Al)sm( 



+ out + p s )- (22) 



The probability amplitude of the process is given by the 
expectation value (Ho), and the Raman intensity of the 
^4^ (Al) mode is given by the square of the coefficient 
of Al {Al) in (H G ) as 



/ G+ K |(^)cos(*) + (^)sin(*)| 2 , 
I G - oc|«>cos(*)-<^) S inW| 2 , 



(23) 



FIG. 2: (a) The strain is at the angle ip a with respect to 
the a;-axis. The direction of strain is parallel to A+. The 
angle between the strain and the incident (scattered) laser 
polarization is denoted by 9i n (Bout)- The axial vector A^ 
corresponds to the G + mode, while Al corresponds to the 
G~ mode, (b) A hexagonal graphene with zigzag edges (top), 
and that with armchair edges (bottom). 



A. Bulk 

Since the polarization of light, e, is parallel to A, the 
electron-photon interaction, A • v, is proportional to e • v. 
Let the polarizations of incident light and scattered light 
be e ln and e out , respectively. Then the effective photon- 
phonon interaction for the phonon mode A q is given by 



H G = e 2 (e out • va z )\A q ■ v q )(e in ■ va z ). 



(19) 



Here, we have omitted to write the electron propagator 
by assuming a resonance Raman process, in which the 



where we have defined VP = 6*i n + 6 ou t + <Ps ^ 

It is important to recognize that only when (v q ) = 0, 
Jq+ and Iq- can have the crystallographic orientation 
dependence of the Raman intensity: 



J G + cx cos 2 (*), 
I G - ocsin 2 (*), 



(24) 



which were used to fit the observed polarization depen- 
dence of the Raman intensity on the crystallographic 
orientation in strained graphene l 7 ! 18 Note that without 
some constraint for v q , the Raman intensity cannot have 
a crystallographic orientation dependence since the elec- 
tronic dispersion is isotropic about the Dirac point. 51 The 
wave function in a periodic graphene does not yield a 
constraint for the axial velocity v q , so that not only (v q ) 
but also (v q ) can take a nonzero value. In this case, we 



cos(9(k) and (v q ) = sin0(k), where 0(k) 



can have 

is the angle between the wave vector k and the fc^-axis. 
After the integral over the variable 6*(k), the intensity 
becomes independent of the angle Vt. Since we have the 
constraint (v q ) = for both the zigzag and armchair 
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orientations, it is naturally expected for the graphene 
sample with rectangle shape enclosed by zigzag and arm- 
chair edges shown in Fig.QJb) that we still have (v%) = 0, 
and that only (u q ) can have a non-vanishing value. Then 
we can reproduce the crystallographic orientation depen- 
dence. As a matter of course, there remains a question 
of whether or not the constraint (v^) — holds in a 
graphene sample with a general edge shape. A further 
discussion on this point will be given in Sec. |V] 

If the Raman intensity of the G band in the bulk of 
graphene does not have a polarization dependence, there 
are in principle two ways to interpret this. One way is to 
assume that the bulk of graphene with edge is identical 
to the "bulk" of a periodic graphene without edge. In 
this case both I G + and Iq- do not have any polarization 
dependence. The other way is to assume that the exis- 
tence of graphene edge gives rise to some constraint for 
the axial velocity in the "bulk" like (v^) = with respect 
to the states participating in the Raman process. In this 
case both Iq+ and I G - do have polarization dependence, 
but the sum of them I G + + I G - does not. Thus, without 
strain, the two kinds of "bulk" can not be distinct. The 
experimental result a 17 i 18 in the bulk of strained graphene 
indicate that the later interpretation is plausible. 

B. Edge 

We now examine the Raman intensity of the G + mode 
and that of G~ mode for the graphene edge. Since we 
have the constraint (t> q ) = for the standing wave, we 
need to modify Eq. (|19| at the graphene edge as 

(H' G ) = e 2 ((e° ut ■ va z y)(A<* ■ ^)(e in ■ va z ). (25) 

In contrast to the effective interaction for the bulk given 
in Eq. (fH?|) , each interaction operator is replaced with 
the expectation value of the operator, by which the in- 
termediate state can satisfy the constraint for the stand- 
ing wave. Physically speaking, this modification assumes 
that the coherence between ingoing and outgoing states 
of the standing wave is strong, so that the intermediate 
state can not transfer into an ingoing state or an outgo- 
ing state independently. This coherence may be weak in 
the bulk, for which case Eq. (fT§|) would become a better 
approximation. We note that regardless of the weakness 
of the coherence in the bulk, the coherence for the initial 
and final states in the Raman process can lead to the 
crystallographic orientation dependences of the G + and 
G - bands. 

Now, with the constraint (v%) — 0, Eq. ([25]) becomes 
(H' G ) = e 2 A^ef(v*)\ (26) 

This is a mathematical expression of the selection rule for 
the G band near the graphene edge in shortened form. 
Since = e in sin((9 in + ip s ), and e° ut = e out sin(0 O ut + 
ip s ) (see Fig. HD, we have with Eq. (dU) that (H' G ) cx 
(A+ cos (fi s — A^_ sin tp„) sin(6*i n + ip s ) sin(6* out +<^ s ). From 



the coefficients of ^4^ and A^_ in this representation, the 
Raman intensity of each mode is given by 

I G + oc cos 2 ((^ s ) sin 2 (6> in + ip s ) sin 2 (0 out -I- tp s ), 

2 2 2 \ / 

I G - cx sin (vJ s )sin (0 in + ip s ) sin (0 ou t + <Ps)- 

It is amusing to note that Iq- (Iq+ ) vanishes when ip s — 
(ips = 90°). Note also that the ratio Iq-/Iq+ depends 
only on the angle ip s , which is in contrast to the case of 
the bulk. 



C. Kohn Anomaly Effect 

Kohn anomaly effect is useful for illuminating the es- 
sential difference between the predictions of the two mod- 
els for the bulk and edge [Eqs. (Q35J) and ([25])]. For the G 
band, the Kohn anomaly effect is caused by the electron- 
hole pair creation from the phonon which is described as 
a vertical transition in the picture of the Dirac cone^r— 
The probability amplitude for the vertical pair creation 
from a phonon mode A q is given by 

M G = ((A«».«V*)t(A*-«V,)) (28) 

This is the formula for the bulk where we assume that the 
"spin" (pseudospin and valleyspin) of the intermediate 
state can be arbitrary (See Fig. [3]). In other words, the 
propagator of electron in the bulk is proportional to the 
identity matrix. In this case, regardless of the character 
of the initial state, we obtain M G = A q -A q from Eq. (j2"5j) . 
This means that both the LO and TO modes in the bulk 
undergo the Kohn anomaly effect. On the other hand, 
the corresponding matrix element for the edge is given 
by 

M^ = <(A q .<AT 2 )t)(A q -uV z }. (29) 

For the case of edge, we assume that the intermediate 
state is given by the standing wave which has the con- 
straint condition for the "spin" . It is easy to show that 
M' G = (A q x (v)) 2 . This leads to the selection rule for 
the G band at edge since we have (v^) = 0. This is the 
cause of the asymmetry that only the LO mode can un- 
dergo the Kohn anomaly effect at both the zigzag and 
armchair edges i 13 ' 14 

IV. D AND 2D BANDS 

In this section we study the polarization dependences 
of the Raman intensities of the D and 2D (G') bands for 
the bulk and edge. It is shown that the Raman intensity 
of the D band in the bulk does not have a polarization 
dependence, while that of the 2D band in the bulk can 
have the polarization dependence, I213 oc (e m • e out ) 2 . It 
is also shown that the Raman intensities of the D and 2D 
bands in strained graphene do not have a crystallographic 
orientation dependence in the bulk. At the edge, these 



6 



Bulk Edge 




Bloch wave Standing wave 



FIG. 3: Kohn anomaly effect at the bulk and edge. In the 
bulk, the intermediate electron- hole state consists of the Bloch 
plane wave, while at the edge, the intermediate state consists 
of the standing wave. The "spin" of the Bloch state can point 
arbitrary direction, while that of the standing wave can not 
have a component parallel to o v t z . 

bands can exhibit the polarization dependence similar to 
that of the G band and also have the crystallographic 
orientation dependence. 

The off-diagonal term in Eq. ([IJ, <pa x (or <p*a x ), rep- 
resents intervalley phonon modes which are responsible 
for the Raman D and 2D bands. Although the D and 
2D bands consist of several phonon modes with different 
wave vectors which depend on the excitation laser energy, 
we examine the Kekule distortion as the representative 
mode. It is straightforward to show that <f> is a constant 
for Kekule distortion, and \4>\ is about three times larger 
than |A q | for the G band£ The latter can explain why 
the intensity of the 2D band is much larger than that of 
the G band. Let us denote = e l9 |0|, then the el-ph 
interaction for the D band Hd is written by 

Hjy = \<t>\a xTe , (30) 

where Tg = t x cos 8 — r y s'm8. Because the matrices a x r x 
and <J X T V do not appear in v and v q , the D band can 
give us new information on the electronic structure that 
is not included in the G band. We will leave the phase 8 
of (f> unspecified because 8 relates to the TO modes near 
the r point through a gauge transformation as is shown 
in Sec. ITVCl 



A. D Band 

The effective photon-phonon interaction for the Kekule 
mode is written by 

H D = H(e out • va z )\a x r e )(^ ■ va z ). (31) 

It is easy to find that 

(H D ) oc e in • e out {a x rg) - ze in x e out (a y r e+f }. (32) 

Hence, when e ln is parallel with e out , i.e., when e m x 
£ out _ ^ e Q band intensity is proportional to 
|(ca;Te)| 2 . On the other hand, when e ln is perpendicu- 
lar to e out , i.e., e in • e out = 0, the D band intensity is 
proportional to \(a y Tg + *)\ 2 . If \{<t x t 6 )\ ± \(a y Tg + *)\, 



the D band intensity in the bulk can have a polariza- 
tion dependence. However, we could not find any special 
reason for this asymmetry. Rather, it is probable that 
\{&xT~e)\ = \{cr y Tg + *)\ holds in the bulk. It is reasonable 
to consider that the D band intensity does not have a 
polarization dependence in the bulk. 

Since the zigzag edge is not the source of intervalley 
scattering, we have (crjT x ) = and (<JiT y ) — for the 
standing wave near the zigzag edge. Thus, we get 

(H D ) = 0. (33) 

This shows that the D band intensity is suppressed near 
the zigzag edge. In contrast, the armchair edge is the 
source of intervalley scattering. In fact, the standing 
wave near the armchair edge is given by 

* k (r) = CV fc ^$ k , (34) 

where C is normalization constant and 4>k is the wave 
function of the pseudospin^iii It is easy to show that this 
wave function reproduces (v x ) = and (v q ) = 0, which is 
consistent with the results obtained in Sec. [II] Moreover, 
pseudospin and valleyspin can be calculated separately 
since the armchair edge preserves the pseudospin. From 
Eq. J35), we get 

(o- x t x ) = (a x ) cos(a - 2k x x), 
{&xTy) = (cr x ) sin(a - 2k x x). 

By using these results, we obtain 

(H D ) = \C\ 2 I \(f>\(a x )cos(8 + a-2k x x)dxdy, (36) 
Js 

for the armchair edge. To summarize, the D band has 
an obvious selection rule; the D band is Raman active 
at the armchair edge while it is not active at the zigzag 
edge. This has been a well-known fact which is useful 
in distinguishing armchair-dominated edge from zigzag- 
dominated edgeiiSr— 

The polarization dependence of the D band intensity 
near the armchair edge is different from that in the bulk. 
In fact, we get from (H' D ) = \<j)\(e out ■ va z )^ (a x Tg){e iu ■ 
v<t z ) that 

(H' u ) cx |0| e ; n 6° ut (^} 3 . (37) 

Let the angle between the armchair edge and the polar- 
ization of the incident (scattered) laser be 0; n (O ou t)- 
Then we can use e^ n = e ln cos 6i n and e° ut = e out cos O out 
in Eq. (|37[) . Thus, the D band intensity at the armchair 
edge behaves according to Id oc cos 2 Oj n cos 2 O ou t which 
is maximum when the polarization of the incident (or 
scattered) light is parallel to the edge. 

Here, let us mention experiments on the polarization 
dependence for the D band intensity at the edge. First, 
the polarization dependence of Id oc cos 2 8j n cos 2 O ou t 
is consistent with the observation for graphite edges by 
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Cancado et a/ j 10 ' 30 and the observations for edges of 
single-layer graphene by You et alH-, Gupta et al^ 
and Casiraghi et alM- Secondly, for edges of single-layer 
graphene, Cong et al~ confirmed that not only the D 
band but also the G band follows Jd,g °c c °s 2 ©in for the 
polarization of the incident laser light, whereas You et 
alr^ and Gupta et al*£ did not observe any polarization 
dependence for the G band. On the other hand, for edges 
of bilayer graphene, Begliarbekov et al. 16 showed that the 
G band intensity had the polarization dependence, while 
the D band intensity did not exhibit any polarization de- 
pendence. The polarization dependence of the D band 
at the armchair edge might be sensitive to the number of 
graphene layers. In fact, Gupta et alM shows that three 
Lorentzian components are necessary to fit the Raman 
spectrum of the D band in bilayer graphene, whereas the 
D band spectra in a single-layer graphene can be well 
fitted by a single Lorentzian component. Theoretically, 
by comparing Eq. (|3"7j) with Eq. (j2f))) . we see that the 
polarization dependence of the D band intensity at the 
armchair edge is identical to that of the G band intensity 
at the armchair edge. We also note that when zigzag and 
armchair edges are randomly distributed along a mixed 
edge, the G band does not show a polarization depen- 
dence^ However, even in this case, the D band should 
have a polarization dependence since there is no coun- 
terpart of the D band which can erase the polarization 
dependence. In the case of the G band, two components 
(LO and TO modes) can coexist in the random edge. 
They have different polarization dependence, so that a 
polarization dependence of the G band may diminish in 
the case of a random edge. 

Equation (|37|) leads to the following polarization de- 
pendence on the crystallographic orientation in strained 
graphene, 

7 D oc sin 2 (6> in + (p s ) sin 2 (6» out + (f s ). (38) 

Thus, from Eq. (gZJ) we find that 7 D oc I G + + I G - . This 
might be one of the most interesting consequence for the 
G and D bands concerning with the armchair edge in 
strained graphene. 



B. 2D Band 

The effective el-ph interaction for the 2D band is given 
by the square of Hd as 

H 2D =Hl = H 2 a T . (39) 

Because the effective interaction is proportional to the 
identity matrix ctqtq, no constraint can affect the el-ph 
matrix element for the 2D band. The effective photon- 
phonon interaction for the 2D band in the bulk is given 
by 

H 2 d = (e° ut • va z )W^\ 2 ovT Q )(e™ ■ va z ). (40) 



Then we have 

(Wan) oc (e hl • e out y T - i(e in x e out )a z r z . (41) 

Note that (ctqtq) = 1 holds for any kind of wave function. 
Furthermore, the condition (<t z t z ) = should be satisfied 
in the absence of a magnetic fields Consequently, the 2D 
band intensity follows I20 (e ln -e out ) 2 in the absence of 
a magnetic field. The Raman intensity of the 2D band in 
bulk is maximum when the incident and scattered polar- 
izations are parallel and minimum when they are orthog- 
onal, which is in good agreement with the experimen- 
tal result by Yoon et al£2- The polarization dependence 
of the 2D band closely resembles that of Rayleigh scat- 
tering 3 - because the effective Hamiltonian for Rayleigh 
scattering is given by 

H K = e 2 (e° ut -v<r z )\e in -va z ). (42) 

Note that is the same as %2D except for the numerical 
factor (coupling constant). In contrast, for the edge, the 
polarization dependence of the 2D band intensity is the 
same as that of the D band: 

<^ D ) oc H 2 e^° u W- (43) 

Note, however, that because the el-ph matrix element for 
the 2D band is given by (-&2d) — \4>\ 2 regardless of the 
orientation of the edge, the 2D band intensity appears 
both at the zigzag and armchair edges. The constraint 
works for the optical transition only, and therefore the 
polarization dependence of the 2D band follows the same 
rule for the G band. 

It is known that Eq. ([3T7)) does not cover the defor- 
mations representing a pentagon or heptagonj 34 ' 35 For 
these topological defects, some combination of a a n be- 
sides a x T x and a x r y can appear. Thus, besides the ap- 
pearance of the Raman peak due to the vibrational dy- 
namics specific to the topological defect^ the D and 2D 
bands can have some information on the presence of the 
topological defect. In fact, it is known that a single pen- 
tagon or a single heptagon gives rise to a mixing between 
K and K' points, leading to a sophisticated topological 
effect on the wave function. Recently, the existence of a 
new type of graphene edge called reczag (reconstructed 
zigzag) has been proposed^ Note that this reczag edge 
consists of a pair of pentagon and heptagon along the 
edge. In this case, the topological effect on the wave func- 
tion is not significant because the topological effect of a 
single pentagon is cancelled by that of a single heptagon. 
In fact, a numerical calculation shows the appearance of 
the edge states near the reczag edge3 This indicates 
that the standing wave near the reczag edge is similar to 
that near the usual zigzag edge. 

C. Gauge Transformation and D' Band 

At first sight, due to the momentum conservation, a 
TO mode with small nonzero momentum (q ^ and 
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q| <C |kp|) may cause an intravalley scattering, but is 
not expected to be relevant to an intervalley scattering. 
Here, in terms of the gauge transformation, we shall show 
that such TO modes do not contribute to intravalley scat- 
tering; rather they can be activated through the interval- 
ley scattering. 

The TO modes with small momentum can be repre- 
sented by the derivative of some scalar function y(r) as 



A q 



Vy>(r), 



(44) 



where A^ on the right-hand side is the zero mode which 
has been relevant to the G band^ Due to the following 
gauge transformation, the scalar function can be trans- 
ferred to the phases of the wave function and <f> as 



<r • (p + A q 



TO J 



o-'-(p-A^o) 



(45) 



where * K = e'^W^K, *K' = e +i ^( r )* K ', and cf> = 
e -2i< ( 3(r)0 j\f t- e t na t t ne XO mode appears as the phase 
of tfi (see 9 in Eq. (f50|) ). A physical significant of this 
gauge transformation is that the TO mode can be excited 
in combination with the intervalley scattering. Since the 
armchair edge enhances the intervalley scattering </>, the 
TO mode with small nonzero momentum <p(r) also can 
be excited near the armchair edge. 

The LO modes with small momentum can be repre- 
sented by the derivative of some scalar function x(r) as 



A2 (r) = Aj + Vx( X (r)e 2 ). 



(46) 



In contrast to the TO mode, A^ (r) can not be gauge 
transformed into the phase of <fi because it has a non- 
vanishing field strength: = V x A£ Q (r) ^ 0. 20 Thus, 
these LO modes are responsible for intravalley scatter- 
ing. The D' band 44 ^ observed slightly above the G band 
in Raman spectra (around 1620cm -1 ) is originated from 
these LO modes. The Hamiltonian for the D' band is 
given by 



H D , = A q (r) • v q . 



(47) 



Note that the "spin" structure is the same as that for 
the G band. Thus, the polarization dependence of the D' 
band follows that of the G band in the bulk. At edge, we 
have (Ht>>) = J s A^(r)^ (r)v^(r)dxdy. The intensity 
behaves as Id' °c \^y Ut ^y\ 2 - The polarization dependence 
of the D' band is the same as that of the G band at edge, 
however, it might be difficult to observe the polarization 
dependence of Id> due to its small intensity. 



V. DISCUSSION AND SUMMARY 

Here, we would like to mention the status of the exper- 
iment on the G band for the graphene edge. Cong et alM- 
conducted a systematic research on edges of single layer 
graphene and found, in particular, that there were two 



orientations of the graphene edge (A-edge and Z-edge) 
which exhibited different behaviors against the polariza- 
tion of the incident laser light. The Raman intensity of 
the A-edge (I a) is enhanced when the polarization be- 
comes parallel to the edge and that of the Z-edge (Iz) is 
enhanced when the polarization becomes perpendicular 
to the edge. By using the angle O between the orien- 
tation of the edge and the polarization of the incident 
laser light, they could fit the observed Raman intensities 
I A (0) and J z (6) with 



J A (6) = a + 6cos 2 6, 
h(&) = c + dsin 2 9, 



(48) 



where a, b, c, and d are fitting parameters. In their ex- 
perimental data, the maximum intensity (a + b or c + d) 
is about two times larger than the minimum intensity (a 
or c), so that a/ (a + b) and c/(c + d) is about 1/2. 

The appearance of these two behaviors for the G band 
is consistent with the selection rule at the graphene 
edge. 14 The A-edge is considered to be armchair dom- 
inant edge and the Z-edge is zigzag dominant edge. The 
ratio of the minimum intensity to the maximum inten- 
sity corresponds to the square of the ratio of the zigzag 
(armchair) part to the armchair (zigzag) part in a mixed 
edge. 14 We thus estimate that the A-edge consists of 60% 
armchair and 40% zigzag edge, while the Z-edge consists 
of 40% armchair and 60% zigzag edge. We can get the 
similar value for the data obtained by Begliarbekov et 
alM- who carried out polarization resolved micro- Raman 
spectroscopy at edges of bilayer graphene. The two ori- 
entations (A-edge and Z-edge) were clearly resolved even 
in bilayer graphene, which also suggests that the num- 
ber of graphene layers does not invalidate the selection 
rule for the G band. Cangado et al££ observed that the 
Raman intensity of the G band for a nanoribbon located 
on top of a highly oriented pyrolytic graphite (HOPG) 
has a strong dependence on the incident light polariza- 
tion. They showed that the Raman intensity is maximum 
when the polarization is parallel to the edge of a nanorib- 
bon. Their result is consistent with the selection rule for 
the armchair edge. A notable point in their experiment 
is that the ratio of the maximum intensity to the min- 
imum intensity was very high. We speculate that the 
nanoribbon located on top of HOPG had rather regular 
armchair edge. 

At this moment, we do not know how to make a clear 
distinction between bulk and edge. In other words, there 
exists no criteria by which we can decide whether the 
Raman process is best described by (He) (Eq. (TH)|) ') or 
(H'q) (Eq. ([25])). Since the electron dynamics in graphene 
is given by massless Dirac equation which is a scale-less 
theory, we consider that bulk of graphene can not be com- 
pletely separated from the edge. To put it in an extreme 
way, there is no bulk region in a nanoribbon with perfect 
regular edge, as well as that there is no edge region in 
a nanotube. The problem is the case of a mixed rough 
edge which might bring a characteristic length scale to 
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TABLE I: Polarization dependences of the Raman intensities for the optical phonon modes in strained graphene. x represents 
absence of a polarization dependence. For "Bulk", #j n (# out ) denotes the angle between the strain and the incident (scattered) 
laser polarization, and <p s is the angle between the direction of strain and the zigzag edge (x-axis). For the "Armchair" and 
"Zigzag" edges, Oi n (0 ou t) denotes the angle between the incident (scattered) laser polarization and the armchair or zigzag 
edge. The polarization dependence of the G band Raman intensity in unstrained graphene is given by I G + + I G - , so that the 
polarization dependence on the crystallographic orientation of strained graphene is lost at the bulk. 







G+ 


G~ 


D 


2D (G') 


Bulk 


cos 2 


(Oin + #out + <fia) 


sin 2 (6» in + 9 out + <p s ) 


X 


(e in ■ e out ) 2 


Armchair 


COS 2 (if s ) 


COS 2 (6 in ) COS 2 (6 out ) 


sm 2 (<p s ) cos 2 (9 in ) cos 2 (6 ou t) 


COS 2 (6in) C0S 2 (6 out ) 


COS 2 (6 in ) C0S 2 (6 O ut) 


Zigzag 


COS 2 '((fis) 


sin 2 (6in) sin 2 (6out) 


sin 2 ^) sin 2 (8i n ) sin 2 (6out) 


X 


sin 2 (9in) sin 2 (6out) 



the scale-less theory. Theoretical estimation of the effec- 
tive length is an important issue should be carried out in 
the near future. Experimentally, according to the Raman 
mapping data by Cong et al.^ the effective region from 
the edge in which the description using the standing wave 
is valid, is about 400 nm, which seems to be comparable 
to the Gaussian laser beam waist ^ 

Let us investigate the property of an eigenstate in the 
interior part of graphene. Since a graphene has the edge, 
the electronic wave function is given by the standing 
wave. Furthermore, the expectation value of the velocity, 
(v), must vanish, so that we have (v x ) — and (v y ) = 
for the standing wave. Note that in the case of nan- 
otubes the velocity around the axis of the tube takes 
nonzero value in general. Suppose that the graphene is 
surrounded by the zigzag edges only [See Fig. &b, top)]. 
Then, (v x ) = and (v y ) = must be satisfied at each 
valley since the zigzag edge is not the source of intervalley 
scattering. Because the velocity and the axial velocity are 
related with each other by v q = vt z , we have (v q ) = 0. 
It is amusing to note that in this special case Eq. ([25)1 
suggests that the resonant G band intensity vanishes. In 
contrast, when graphene is surrounded by the armchair 
edges only [See Fig. [2jb, bottom)], we have ^ and 
(v q ) = since the armchair edge is relevant (irrelevant) to 
the valleyspin (pseudospin) . Note that in both hexagonal 
graphenes, we have at least the condition (u q ) — 0. It is 
our speculation based on the observation of several edge 
shapes that the constraint (v q ) = holds for graphene 
with a wider variety of edge shapes. 

The formalism using the gauge fields for photon A and 
phonon A q might provide a new and fresh insight into 
Raman scattering, which otherwise well-studied subject. 
In ordinary Raman spectroscopy, we irradiate a laser 
light A m onto a graphene sample and observe the in- 
elastically scattered light A out . This Raman process in 



unstrained graphene may be represented as 

A in^ A q + A out_ ( 49 ) 

The left-hand side of "— shows the input and the right- 
hand side of it denotes the output. Note that a phonon 
A q on the right-hand side is a kind of lattice deforma- 
tion or an internal strain. Thus, in Raman spectroscopy, 
by inputting a photon (electronic) signal, one gets a sig- 
nal of strain from graphene. Let us consider a process 
represented by 

A? n ^A + A q ut , (50) 

where A q represents an external strain. This process 
of Eq. ([50)1 is given by replacing A (A q ) with A q (A) in 
Eq. (|49[) . In this process, by inputting strain, one gets an 
electronic output A from graphene, which seems to be a 
prototypical process of strain engineering. Now, the Ra- 
man process in strained graphen o 17 i 18 ' 41 ~ — is expressed 
by 

AL k + A in ^AS ulk + A q + A° ut , (51) 

where the strain is described by A^ ulk . Considering 
that this is a process which may be recognized as the 
sum of Eqs. (|4"9")l and ([5D|) . it can be said that Raman 
spectroscopy in strained graphene is a small step toward 
strain engineering in graphene. 

It is reasonable to consider that Eq. (|49]) represents a 
Raman process in graphene without edge. Strictly speak- 
ing, a Raman process in a real (unstrained) graphene 
should be represented not by Eq. (l49l) but by 

A? dgc + A i "^A^ dgc + A q + A- t , (52) 

because the presence of the edge is represented by a 
local strain field A^ dgo : the zigzag edge is represented 
by a local A q field parallel to the edge, while the arm- 
chair edge corresponds to a local A q field normal to the 
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edge^^ It is interesting to note that the direction of 
the A*L. is coincident with that of the Raman active 

CQ & C 

phonon mode A q near the graphene edge. Considering 
that Ag ulk in Eq. d5j) (A|? dgc in Eq. ([52])) represents a 
global (local) strain, the Raman spectroscopy near the 
graphene edge^2 r 12 i 15 ' 16 is complementary to the Raman 
spectroscopy in the bulk of strained graphene. 

In conclusion, the polarization dependences of the G, 
D, and 2D Raman bands at bulk and edge have been 
investigated theoretically with paying attention to the 
pseudospin and valleyspin of the standing wave. Our re- 
sults are summarized in TABLE [J] The constraint for 
the axial velocity provided by the graphene edge is es- 
sential to the selection rule for each Raman band. The 
selection rules of the G and D bands for graphene edge 
are consistent with the recent experimental results. The 
coherence provided by the graphene edge seems to persist 
even in the bulk, by which we explain the recent experi- 
ments for strained graphene showing the crystallographic 



orientation dependences of the Raman intensities of the 
G + and G~ bands. This also suggests that the "bulk" of 
graphene can not be completely free from the surround- 
ing graphene edge, and that it is necessary to distinguish 
the bulk of graphene that is surrounded by edge from the 
bulk of a periodic graphene without edge. 
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tial in the subsequent arguments. 

50 In literature, 3<p s instead of <p s appears in the angle 'f. The 
factor of 3 in front of ip„ seems to come from the choice of 
the coordinate system, for which tp s in Eq. (|21|) is replaced 
with — </5 s .— We could not find any proper reason for this 
choice. 

51 This feature of graphene bulk is different from that of nan- 
otube because the electronic dispersion is not isotropic in 
the case of nanotube due to the cutting lines. 



